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A fundamental gap in the study of the origin of limbed vertebrates lies in 
understanding the morphological and functional diversity of their closest relatives. 
Whereas analyses of the elpistostegalians Panderichthys rhombolepis, Tiktaalik roseae 
and Elpistostege watsonihave revealed a sequence of changes in locomotor, feeding 


and respiratory structures during the transition! °, anisolated bone, a putative 
humerus, has controversially hinted at a wider range in form and function than now 
recognized’”* “. Here we report the discovery of anew elpistostegalian from the Late 
Devonian period of the Canadian Arctic that shows surprising disparity in the group. 
The specimen includes partial upper and lower jaws, pharyngeal elements, a pectoral 
fin and scalation. This new genus is phylogenetically proximate to 7. roseae and 

E. watsoni but evinces notable differences from both taxa and, indeed, other 
described tetrapodomorphs. Lacking processes, joint orientations and muscle scars 
indicative of appendage-based support ona hard substrate”, its pectoral fin shows 
specializations for swimming that are unlike those known from other sarcopterygians. 
This unexpected morphological and functional diversity represents a previously 
hidden ecological expansion, asecondary return to open water, near the origin of 


limbed vertebrates. 


Study of tetrapodomorph skulls, fins, axial skeleton and scalation 
has revealed the ways that feeding, respiration and appendage-based 
locomotion changed as fish shifted from aquatic to terrestrial life- 
styles>"*. P. rhombolepis*®, T. roseae*® and E. watsoni’ hold a special 
place in these analyses, showing a combination of plesiomorphic and 
apomorphic features that give insight into a sequence of anatomical 
changes in the origin of limbed taxa (that is, those in possession of 
digited appendages and lacking dermal rays). Now missing, however, 
is an understanding of the morphological, functional and ontogenetic 
diversity of the finned tetrapodomorphs most closely related to limbed 
forms. This is unfortunate, as isolated or fragmental specimens have 
controversially hinted at a wider range of diversity thanis observed in 
more complete material”. 

Here we describe a new finned tetrapodomorph thatis closely related 
to T. roseae and E. watsoni. The new form shows an unexpected com- 
bination of characters, one that suggests a broad range in disparity 
among the closest finned relatives of limbed forms. The specimen was 
collected 1.5 km east of the site that yielded T. roseae, but fromaslightly 
lower horizon in the Fram Formation of southern Ellesmere Island, 
Nunavut, Canada. We describe this new taxon and present a phyloge- 
netic analysis to reveal its implications for understanding the evolution 
of the nearest relatives of limbed tetrapodomorphs. Comparison of 
the new taxon to other Frasnian-age forms allows a reinterpretation 
of isolated elements of previously uncertain affinity, thus, indicating 


amore widespread and diverse assemblage of tetrapod relatives than 
previously recognized. 


Geological framework 


Embry and Klovan” described the type section of the Fram Forma- 
tion from a drainage feeding the eastern arm of Bird Fiord on south- 
ern Ellesmere Island. They indicate an early to middle Frasnian age 
for the Fram Formation on the basis of palynological spot samples, 
which were collected from near the base, the middle and top of the 
formation”. The Nunavut Paleontological Expeditions collected ver- 
tebrate remains from 2000 to 2008 at 16 sites from the Fram Forma- 
tion within the type section. The holotype of T. roseae (NUFV 108), as 
wellas all other 7. roseae specimens, were collected from site NV2K17, 
which occurs in silty overbank floodplain deposits'® at 533 m above 
the base of the measured type section of Embry and Klovan”. The 
specimen discussed here (NUFV 137) was collected at site NVO401 
(77° 10.235’ N, 86° 11.279’ W) from lower in the same section and 
1.5 km from NV2K17 (Fig. 1a,b and Extended Data Fig. 1). Site NVO401 
is about 453 m above the base of the type section and occurs ina 
medium-grained sandstone. The surface-collected NUFV 137 is the 
only specimen found at the site. NUFV 137 is older than T. roseae and 
was collected from a different facies in the floodplain deposits of the 
Fram Formation. 
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Fig. 1|Locality and holotype of Qikiqtania wakei gen. et sp. nov. a, Specimen 
NUFV 137 was discovered on southern Ellesmere Island, Nunavut, Canada. 

b, The site, NVO401, lies 80 m below NV2K17, the site where T. roseae was 
discovered. c, Materials were pCT scanned and are shown here in dorsal aspect. 
General body shape based on specimen MHNM 06-2067 of E. watsoni’. Scale 
bar, 10 cm. Abbreviations: Fm, formation. 


Systematic palaeontology 


Sarcopterygii’” Romer, 1955 
Tetrapodomorpha” Ahlberg, 1991 
Elpistostegalia™ Camp and Allison, 1961 


a 
sbp 


Fig. 2 | The feeding apparatus of Q. wakei. Volume renderings of CT scans 

of the lower jaw and extra fragments reconstructed in their natural positions. 
a, Dorsal view of the lower jaws, ceratohyal, gular plates, premaxilla and palate. 
Scale bar, 1 cm. b, Ventral view with premaxillary and palatal elements 
displaced so ventral surfaces are visible. c, Left lower jaw, dorsal. d, Left lower 
jaw, medial. e, Right lower jaw, ventral. f, Right lower jaw, lateral. Abbreviations: 
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Qikiqtania wakei gen. et sp. nov. 

Locality. Canada, Nunavut, southern Ellesmere Island, near the 
eastern arm of Bird Fiord, Nunavut Paleontological Expedition site 
NV0401, 77° 10.235’ N, 86° 11.279’ W. 

Geological setting. Fram Formation (Upper Devonian, early Frasnian 
Stage). 

Etymology. Qikigtania (pronounced ‘kick-kiq-tani-ahh’) is derived 
from Inuktitut word Qikiqtaaluk/Qikiqtani, the traditional name for 
the region where the fossil site occurs. The species designation is in 
memory of David Wake, an eminent evolutionary biologist and trans- 
formative mentor, late of the University of California at Berkeley. 

Holotype. Nunavut Fossil Vertebrate Collection (NUFV) 137. 

Material. The description is based ona specimen from the NV0401 
site that preserves the lower jaws, partial left upper jaw and palate in 
articulation, gulars, ceratohyals, an articulated left pectoral fin and 
articulated scales from the dorsal midline, flank and lateral line series 
(Fig. Ic, Extended Data Fig. 2 and Supplementary Video 1). The jaw 
material was physically prepared at the Academy of Natural Sciences 
of Drexel University. Computed tomography scans were collected at 
The University of Chicago’s PaleoCT scanning facility (Supplemen- 
tary Table 1). The specimen will be housed at the Canadian Museum 
of Nature, Ottawa, Ontario, until such time as research and collections 
facilities are available in Nunavut. 

Diagnosis. Elpistostegalian tetrapodomorph characterized by the 
following unique combination of characters: dorsoventral asymmetry 
in pectoral fin lepidotrichia (also present in 7. roseae) and possession 
of aboomerang-shaped humerus lacking ventral ridge and associated 
foramina and ectepicondyle (distinct from P. rhombolepis, F. watsoni, 
T. roseae and more crownward tetrapods). 


Description 

Upper jaw and palate 

Rostral elements of the upper jaws and palate, including portions of 
the ectopterygoid, dermopalatine, vomer, premaxilla and maxilla are 
preserved (Fig. 2a,b, Extended Data Fig. 2and Supplementary Video 2). 
These elements are primarily from the left side and preserved in articu- 
lation with the lower jaws. The vomer is broad, fanged and forms the 
posterior wall of the palatal fossa with a row of smaller teeth. Fangs 
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ac, anterior coronoid; acf, anterior coronoid fang; ch, ceratohyal; cho, choana; 

d, dentary; df, dentary fang; dpf, dermopalatine fang; ecf, ectopterygoid fang; 

g, gulars; mc, Meckel’s cartilage; mcf, Meckelian canal foramen; mx, maxilla; pa, 

prearticular; pc, posterior coronoid; pcf, precoronoid fossa; pmx, premaxilla; 


pspl, postsplenial; sbp, submandibulo-branchiostegal plate; psf, 
post-symphyseal flange; vf, vomerine fang. 


Fig. 3 | Left pectoral fin of Q. wakei. Volume renderings of CT scans of the 
fin with scales removed. a,b, Dorsal (a) and ventral (b) views of the fin with 
endoskeleton in grey and dermal rays in orange. Dotted lines indicate the 
boundary between ulna and ulnare. The dashed line indicates position of 
cross-section ine, which is oriented orthogonal to the plane of the fin web. 
c, Endoskeleton viewed from the proximal side with humerus removed. 

d, Reconstruction of endoskeletal elements with estimated boundary 


and a row of smaller teeth are also present on the dermopalatine and 
ectopterygoid. An expanded mesial surface of the dermopalatine lacks 
teeth and overlaps slightly with the vomer, similar to 7. roseae®, form- 
ing the mesial and posterior margin of the choana. The anterolateral 
wall of the choana is formed by a simple, smooth articulation of the 
premaxilla and maxilla. Maxillary teeth are smaller than the premaxil- 
lary teeth. In their respective tooth rows, maxillary and premaxillary 
teeth are uniform in size. 


Lower jaw 

The lower jaws of Q. wakeiare preserved in articulation anterior to the 
adductor chamber, including the dentary, infradentaries, coronoids 
and prearticular (Fig. 2). The symphysis is relatively smooth, not inter- 
digitating. Large fangs with plicidentine infolding are present on the 
dentary, anterior coronoid and middle coronoid. Rows of smaller denti- 
tion are also present onthe coronoids and dentary, including evidence 
of an auxiliary lateral tooth row on the dentary. The prearticular has a 
broad shagreen field of denticles that is raised adjacent to coronoids, 
and the denticles possess a distinct dorsoventral gradient in size. The 
adsymphyseal is missing, but small teeth embedded in the matrix of 
the precoronoid fossa suggest it was present in life. 


between the radius and intermedium. Scale bar, 1cm.e, Cross sections of 

the fin rays, showing asymmetry in the size of dorsal and ventral hemitrichia. 
f-k, Humerus in dorsal (f), pre-axial (anterior) (g), ventral (h), postaxial 
(posterior) (i), proximal (j) and distal (k) perspectives. Scale bar, 1 cm. Proximal 
is up inf-i. Dorsal is up injandk. Abbreviations: ar, anterior radial; cap, caput 
humeri;h, humerus; ir, intermedium; r, radius; rf, radial facet; m3, third 
mesomere, u, ulna; ul, ulnare; uf, ulnar facet. 


Infradentaries are identifiable by the presence of the mandibular 
canal and postsplenial pit line. The mandibular canal is an open groove 
along most of its length, but in areas of the most intact preservation 
it takes the form of discrete pits the bone surface. The splenial has a 
larger postsymphyseal flange than in T. roseae but has a similar articula- 
tion with the prearticular*. Boundaries between the infradentaries are 
obscured by overlying dermal sculpting and are difficult to distinguish 
in computed tomography cross-section. 

The Meckelian canal contains only partially ossified Meckelian bone 
along its length, but evidence of Meckelian ossification extends from 
the symphysis to the posterior coronoids. The canal is exposed lingually 
ventral to the prearticular and, in areas of intact ossification, Meckelian 
fenestra are bordered dorsally by Meckelian bone and ventrally by 
infradentaries. 


Gular plates and ceratohyal 
Fragments of a principal and median gular plate are preserved, along 
with a series of submandibulo-branchiostegal plates (Fig. 2a,b). A 
straight, grooved ceratohyal lies immediately adjacent to the left lower 
jaw. Asmall anterior fragment of the right ceratohyal is preserved adja- 
cent to the right lower jaw? 7. 
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Fig. 4| Phylogenetic analysis. a, Strict consensus tree from the maximum 
parsimony analysis with Bremer decay (D) and bootstrap support values. 
b, Majority-rule tree from undated Bayesian analysis with posterior 
probabilities. Both analyses recover a basal polytomy; Megalichthys is 
shownas the outgroup, consistent with other studies®”*?>*, 


Pectoral fin 

The left pectoral fin includes the humerus, ulna, radius, intermedium, 
third mesomere, third radial, fin web and associated scales (Fig. 3a,b 
and Supplementary Video 3). The fin is embedded in matrix with the 
proximal articular surface of the humerus and the posterior distal fringe 
of the fin web exposed at the edges of the block (Extended Data Fig. 3a). 
Three endoskeletal elements contact the humerus. Two have robust 
proximal articular surfaces and are identified as the radius and ulna. 
The third, which lies between and slightly dorsal to them, is identified as 
the intermedium proximally displaced during preservation, although 
its shape is difficult to assess because of its position relative to other 
elements (Fig. 3c,d and Supplementary Discussion). 

The fin is characterized by ventralward curvature of the radius and 
asymmetry in the lepidotrichia, in which dorsal hemitrichia have a 
greater cross sectional area than ventral hemitrichia, as in T. roseae 
(Fig. 3e and Extended Data Fig. 3c,d)’. Roughly 30 lepidotrichia are 


Dorsal 


Ventral 


Tiktaalik 


Eusthenopteron 


Panderichthys Qikiqtania 


preserved. Similar to other finned tetrapodomorphs, rays are more 
robust anteriorly and more gracile posteriorly and rays are more ter- 
minally positioned onthe posterior side’. 

The humerus is boomerang-shaped and lacks numerous charac- 
teristic elpistostegalian features, notably a humeral ridge and asso- 
ciated foramina, ectepicondylar process, prominent entepicondyle 
and distinct articular surfaces for the ulna and radius (Fig. 3f-k). The 
ulna lacks a postaxial process and distally would have articulated with 
the intermedium and ulnare. The fin is gracile as compared to other 
elpistostegalians. The anteroposterior width of the humerus is nar- 
rower than the humeriof T. roseae’ and £. watsoni? and more similar to 
P. rhombolepis’. The shallow dorsoventral depth of the fin might reflect 
compression; however, articular surfaces of the ulna and radius are 
similar in their geometry to three-dimensionally preserved specimens 
of T. roseae specimens (NUFV 108, 109, 110), suggesting that morphol- 
ogy was narrow in life (Supplementary Discussion). 


Scalation 

Scales are preserved from the trunk, including dorsal midline and flank, 
the pectoral fin and the lateral line series (Extended Data Fig. 4). Scala- 
tionis broadly similar to other finned elpistostegalians”’””. Scales are 
rhomboid in shape with the free surface sculpted and a smooth internal 
surface that often bears a ventral keel (Extended Data Fig. 4a—c). Onthe 
trunk, scale rows extend posterolaterally from the dorsal midline, with 
individual scales partially covering the scale that follows inthe row and 
also the scale of an adjacent posterior row (Extended Data Fig. 1d,e). 
Pectoral fin scales are smaller than those of the flank and show varia- 
tion in their morphology (Extended Data Fig. 4f-m). Lateral line scales 
are preserved from the left flank and showa completely enclosed tube 
with anterior suprascalar and posterior infrascalar pores enlarged 
relative to the diameter of the canal, and a small pore midway along the 
length of the scale connecting the canal to the external environment 
(Extended Data Fig. 4n-r). 


Phylogenetic relationships 

The phylogenetic position of Q. wakei was analysed by maximum par- 
simony and undated Bayesian approaches, which were applied to a 
matrix of 13 taxa and 125 characters primarily assembled from previous 
publications (see Supplementary Information for detailed description 
of phylogenetic data)?”*. Both methods robustly recover Q. wakei as 
crownward to P. rhombolepis and, thus, as an elpistostegalian closely 
related to limbed tetrapods (Fig. 4). The analyses differ in their relative 


Elpistostege Acanthostega GSM 104536 


|________| 


Ichthyostega 


Fig. 5| Humeriat the fin-to-limb transition. For consistency of orientation 
between species, several specimens have been reflected, so that each is 
represented as being from the right side. Illustrations are based on previously 
published descriptions: Fusthenopteron”®, Panderichthys*”*, Tiktaalik°, 
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Elpistostege’, Acanthostega’’, Ichthyostega®® and GSM104536 (refs. *°"). 
Abbreviations: ect, ectepicondyle; ent, entepicondyle; hr, humeral or ventral, 
ridge. 


placement of Q. wakei, T. roseae, F. watsoni; a strict consensus tree of 
the 28 shortest trees recovered from maximum parsimony analyses 
shows an unresolved polytomy, whereas Bayesian analysis finds weak 
support for a sister relationship between Q. wakeiand T. roseae with 
E. watsoni positioned more crownward. This is similar to other recent 
phylogenetic analyses of stem tetrapods, which have robustly recov- 
ered Tiktaalik and Elpistostege as outgroups to digited forms, although 
support for their relative positions is not strong?”*”>. 


Discussion 


Q. wakei reveals a combination of characters unique among stem tet- 
rapods. The pectoral fin, lacking a postaxial process onthe ulnare and 
showing accentuated hemitrichial asymmetry, is clearly elpistostega- 
lian®”. Yet, the morphology of the humerus is unlike others described. 
With the absence of a ventral ridge or ectepicondylar process and in 
possession of a general boomerang shape, it is more similar to the 
humerus previously attributed to the tetrapod, Flginerpeton pancheni, 
than to any other Devonian taxon (Fig. 5). That specimen, GSM 104536, 
from Scat Craig in Scotland, is an isolated bone from a coeval deposit 
in Laurentia that generated debate as to whether it was from atetrapod 
or whether it was even a humerus at all’? “. The similarity to Q. wakei 
suggests that GSM 104536 is indeed a humerus but belongs toa finned 
elpistostegalian, not a limbed tetrapod. 

The morphology of the Q. wakei humerus is distinctive among stem 
tetrapods. Indeed, the lack of muscular processes on the humerus for 
flexors and extensors at the shoulder and elbow, the terminal position 
of the facets for the radius and ulna, and the relatively large surface 
area of the fin web suggest that the fin of Q. wakei is less suited for 
walking, trunk lifting and station holding in water than itis for a range 
of swimming behaviours”. With its gracile form and lacking many of 
the known main osteological correlates of muscular attachment”®, the 
pectoral fin of Q. wakei represents a strategy of controlling hydrody- 
namic forces not seen in other stem tetrapods. As these features are 
not seen in tristichopterids, osteolepids or rhizodontids, they probably 
arose as apomorphies in elpistostegalians. 

The holotype of Q. wakei is estimated to be a standard length of 
75 cm (calculated from the proportions of £. watsoni specimen MHNM 
06-2067 (ref. °) scaled to the length of the lower jaw), making it smaller 
than other described elpistostegalians. The ontogenies of Eusthenop- 
teron foordiand T. roseae provide evidence that, despite its relatively 
small size, the unique humeral morphology of Q. wakei reflects phylo- 
genetic signal and not developmental stage. F. foordiindividuals are 
described spanning more than 40-fold variation in size”’, and acrossa 
broad range of sizes uniformly retain a ventral ridge, entepicondylar 
process and orientations of facets for articulation with the radius 
and ulna”®”’. 7. roseae, which is known from humeri ranging twofold 
in size, show a similar pattern, preserving these features across this 
size range, although overall proportions might vary*’”. Thus, major 
ontogenetic shifts in limb skeletal anatomy of /chthyostega and Acan- 
thostega, indicated to correspond to aquatic subadults transitioning 
to more terrestrial adult lifestyles using appendage-based substrate 
support, are derived for limbed forms*°. Finned tetrapodomorphs, 
by contrast, are predicted to show more minor changes in the pro- 
portions of endoskeletal, and potentially dermal, components of 
their paired fins’. 

With two elpistostegalian genera now known from nearby localities 
in Canadian Arctic and others from Quebec’, Latvia”’” and poten- 
tially Russia®’, Australia** and Scotland”, the group probably has a 
wide distribution by the Frasnian Stage of the Late Devonian. This 
broad biogeographic range, coupled with the morphological disparity 
revealed by Q. wakei, hints at a wider diversity of elpistostegalians 
than knownat present, with the closest relatives of tetrapods adapt- 
ing in new ways to benthic, littoral and open water habitats by the 
Late Devonian”>*, 
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Methods 


Computed tomography scanning 

LCT scans were collected at The University of Chicago’s PaleoCT 
scanning facility using a GE Phoenix v|tome|x 240 kV/180 kV scanner 
(http://luo-lab.uchicago.edu/paleoCT.html). Scan parameters are 
reported in Supplementary Table 1. w.CT data were reconstructed 
with Phoenix Datos|x 2 (v.2.3.3), imported to VGStudio Max (v.2.2) 
for cropping and exportation as a 16-bit tiff stack. Tiff stacks were 
segmented and visualized in Amira v.20.2 (FEI Software). For some 
scans, to accommodate for computational challenges that arise from 
large file sizes, data were converted to 8-bit files for segmentation; in 
such cases, after segmentation the renderings were generated from 
the original 16-bit files. Animations were generated by exportation 
tiff stacks from Amira and then edited with Adobe Premiere (v.13.12). 
High-resolution versions of images from all figures are provided in 
the Supplementary Data. 


Phylogenetic analyses 

We investigated the phylogenetic position of Q. wakei using a phy- 
logenetic data set of 13 taxa and 125 characters (see Supplementary 
Information for a detailed description of the phylogenetic matrix). 
All characters were treated as equally informative, and we assumed 
unordered evolution among states. 

Maximum parsimony analyses were performed using PAUP* 
(v.4.0a168)**. The branch and bound method for searching tree 
space was used with no topological constraints. A total of 28 
most-parsimonious trees were recovered (tree length 151s). The trees 
are summarized as a strict consensus tree (Fig. 4) and as an Adams con- 
sensus tree (Extended Data Fig. 5a). Clade support was estimated using 
two approaches: Bremer decay values”, calculated with AutoDecay 
(v.5.06)*° and non-parametric bootstrapping, calculated in PAUP* with 
500 replicates (Fig. 4 and Extended Data Fig. 5b). Apomorphies of nodes 
inthe strict consensus tree were identified using the function ‘apolist’ 
in PAUP*, which returns characters optimized under both accelerated 
transformation (ACCTRAN) and delayed transformation (DELTRAN) 
conditions (Extended Data Fig. 5c). 

Undated Bayesian analyses were performed using MrBayes 
(v.3.2.7a)*. Analyses were run for five million generations with four 
runs of four chains sampling every 5,000 generations and a burn-in 
of 20%. Megalichthys was designated as an outgroup, consistent with 
other studies’”*”5*°, 

Convergence was assessed with diagnostics reported by MrBayes 
(average s.d. of split frequencies <0.02, potential scale reduction fac- 
tors was 1, effective sample sizes >200). Results are summarized by a 
majority-rule consensus tree of postburn-in trees (Fig. 4). 

For both maximum parsimony and Bayesian analyses, executable 
files, log files, and individual trees that contribute to the summary trees 
are included as supplementary files (Supplementary Data). Custom 
R (v.3.6.1) code used for the calculation of Bremer decay values and 


for visualization of phylogenies are available at https://github.com/ 
ThomasAStewart/Qikiqtania. All code is archived at Zenodo (https:// 
doi.org/10.5281/zenodo0.6557684). 


Reporting summary 
Further information on research design is available in the Nature 
Research Reporting Summary linked to this paper. 


Data availability 


All data are available for download. Computed tomography data sets 
and STL files of main elements are available for download from Mor- 
phoSource (https://www.morphosource.org/projects/000375542). 
Phylogenetic data are in Supplementary Data Files 2 and 3. Source data 
are provided with this paper. 


Code availability 


Executable files for PAUP* and MrBayes are available in the supple- 
mentary materials. Custom R code used for the calculation of Bremer 
decay values and for visualization of phylogenies are available at 
https://github.com/ThomasAStewart/Qikiqtania. All code is archived 
at Zenodo (https://doi.org/10.5281/zenodo.6557684). 
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Extended Data Fig. 1| Photograph of the locality NVO401. Photograph arrows highlight approximate stratigraphic separation between the two 
showing the localities NVO401, where NUFV 137 was collected, and NV2K17, horizons. White lines trace two additional horizons across the valley. A yellow 
where 7. roseae was collected. White arrows indicate sites of collection. Yellow tent, approximately 2.5 macross, isinthe midground. 


Extended Data Fig. 2 | Photographs of NUFV 137. (a-f) Elements associated 
with the feeding apparatus. Elements in a-d are shown in Fig. 2 and Video S2. 
Element eis identified as parts of the palate and lower jaw due to the presence 
of multiple rows of both dorsally and ventrally facing teeth. The ventrally facing 
teeth are determined to be palatal in nature due to the expanded medial 
shagreen of denticles (likely part of the entopterygoid) that are bordered 
laterally by two uniform rows of larger teeth (likely the ectopterygoid and 
maxilla). This piece could not be definitively positioned relative to the other 


jaw elements due to absence of the corresponding broken tooth bases onthe 
main lower jaw block. Element fis identified as part of a lower jaw on the basis of 
its curvature and dentition. (g) Maxilla. (h) Left pectoral fin, whichis 
embedded in matrix, with exposed associated scales. (i) Fragment containing 
scales and lepidotrichia froma paired fin. (j, k) Fragments with undiagnosed 
vascularized endoskeletal elements. (I) Scale field from the dorsal midline, 
anterior is up. (m) Fragment containing scales from the lateral line series and 
flank. (n) Trunk scale field, anterior is left. 
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Extended Data Fig. 3 | Additional fin-associated materials. (a) Athin, slightly —lepidotrichia froma paired fin. (c) Finrays, seen in the lower right corner of 
convex bladelike element that might be part of the pectoral girdleisadjacentto panelb,intheir preserved position showing asymmetry between the dorsal 
the pectoral fin. Breaks in the block have exposed the proximal articular and ventral hemitrichia. (d) Three pairs of hemitrichia from panel c 

surface of the humerus and the posterodistal portion of the fin web. (b) An repositioned and shown in dorsal perspective. Dorsal hemitrichia are orange, 
element, shown in Extended Data Fig. 2i, contains scales and additional and ventral hemitrichia are blue. 


Extended Data Fig. 4| NUFV 137 scales. Internal (a) and external (b) views of 
scale field from left flank. (c) Scales outlined in panel b showing median ridge 
oninternal surface. Internal (d) and external (e) views of scale field from dorsal 
midline. (f) Left pectoral fin in ventral aspect showing the position of individual 
scales figured in panels g-l. (g) Scales that covered the humerus ventrally. (h, i) 
Elongate scales from leading edge. (j-I) Small scales from the ventral surface of 
the fin. (m) One scale in pre-axial (anterior), external, post-axial (posterior) and 
internal views showing dermal sculpting and lack of ventral keel. (n, 0) Left 


lateral line scale in external and internal views. (p, q) Scale with reduced opacity 
and the canal shown in blue. Midway along the length of the scale, a pore 
connects fluid in the canal and the external environment. (r) Two of the 
preserved lateral line scales in reconstructed position showing their degree of 
overlap and expected orientation relative to the epidermis. In panels m,n: area 
of overlap with adjacent scale in the row shown in orange, area of overlap with 
scale in adjacent row shown in pink. Abbreviations: asp, anterior suprascalar 
pore; pip, posterior infrascalar pore; p, pore. 
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a Maximum parsimony, b Maximum parsimony, 
Adams consensus Majority-rule of 500 bootstrap 

Megalichthys Megalichthys 
Tinirau Tinirau 
Eusthenopteron Eusthenopteron 
Panderichthys ad Panderichthys 
Qikiqtania a Tiktaalik 
Tiktaalik x) 7 Qikiqtania 
Elpistostege i Elpistostege 
Ymeria ot Parmastega 
Parmastega -— Ventastega 

— Ichthyostega i 39L_ Acanthostega 

L_ Elginerpeton 28| -——— Elginerpeton 

— Ventastega 21 -— Ymeria 

L_ Acanthostega ~ 39 Ichthyostega 


c Maxi Fr 1 Char. 4 (Frontal): 0 to 1 
aximum parsimony, , Char. 43 (Angular-prearticular contact):1 to 0 
Strict consensus and unambiguous character changes Char. 83 (Intracranial joint): 0 to 1 
Char. 103 (Ribs, length): 0 to 1 
Char. 105 (Ribs, shape): 0 to 1 


Megalichthys Char. 107 (Scapulocoracoid): 0 to 1 
Tini Char. 111 (Postaxial process on ulnare): 0 to 1 
inirau Char. 115 (Long basal segments of lepidotrichia): 0 to 1 
Char. 116 (Basal scutes on fins): 1 to 0 
Eusthenopteron Char. 117 (Tooth construction): 0 to 1 
Panderichthys 2 Char. 2 (Ectopterygoid/palatine exposure): 0 to 1 
Qiki tani Char. 12 (Median rostral): 0 to 1 
ikigtania Char. 14 (Prefrontal, length): 0 to 1 
o ij Char. 15 (Prefrontal, suture): 0 to 1 
Tiktaalik Char. 22 (Basipterygoid process): 0 to 1 
7 Char. 23 (Ethmoid ossification): 0 to 1 
Elpistostege Char. 62 (Prearticular suture): 0 to 1 
Parmastega Char. 63 (Prearticular ridge): 0 to 1 
9g Char. 101 (Posttemporal + supracleithrum): 0 to 1 
Ventastega Char. 110 (Proportion of skull roof): 0 to 1 
i, 3 Char. 8 (Lateral rostral present: 0 to 1 
Acanthostega Char. 17 (Pterygoids separate): 0 to 1 
3 ‘ Char. 25 (Otic capsule): 0 to 1 
Elginerpeton Char. 44 (Coronoid (anterior) contacts splenial): 0 to 1 
Ymeria Char. 47 (Dentary external to angular + surangular): 0 to 1 
Char. 52 (Meckelian bone floors precoronoid fossa): 0 to 1 
L__ Ichthyostega Char. 57 (Adsymphysial mesial foramen): 0 to 1 


Char. 69 (Coronoid fangs): 0 to 1 

Char. 73 (Dentary teeth: larger than maxillary teeth): 1 to 0 
Char. 74 (Dentary with a row of very small teeth): 0 to 1 
Char. 76 (Adsymphysial plate dentition): 0 to 1 

Char. 77 (Adsymphsial plate fang pair): 0 to 1 

Char. 78 (Adsymphysial plate tooth row): 0 to 1 

Char. 79 (Prearticular shagreen field): 0 to 1 

Char. 84 (Nature of dermal ornament): 0 to 1 

Char. 85 (Nature of ornament): 0 to 1 

Char. 90 (Digits): 0 to 1 

Char. 91 (Humerus, entepicondy): 0 to 1 

Char. 92 (Pectoral process of humerus): 0 to 1 

Char. 100 (Lepidotrichia in paired appendages): 0 to 1 
Char. 107 (Scapulocoracoid): 1 to 2 

Char. 113 (Sacrum): 0 to 1 

Char. 118 (Gular): 0 to 1 


Extended Data Fig. 5| Expanded results of phylogenetic analyses. (a)Adams_ _ phylogenetic analyses of early tetrapods?”*”>*°, although basal polytomies are 
consensus tree of maximum parsimony analyses. (b) Majority rule tree of recovered. (c) Unambiguous character changes recovered onthestrict 


maximum parsimony analyses with bootstrapping (500 replicates). In all consensus tree using the command ‘apolist’ from PAUP**®. 
panels, Megalichthys is plotted as the outgroup consistent with previous 
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Policy information about availability of computer code 


Data collection CT data were collected on a GE Phoenix v|tome|x 240 kv/180 kv scanner, reconstructed with Phoenix Datos|x 2 (v2.3.3), and then imported 
to VGStudio Max (v2.2) for cropping and exportation as a 16-bit tiff stack. 


Data analysis Tiff stacks were segmented and visualized in Amira v20.2 and Adobe Premiere (v13.12). Code for PAUP*(v4.0a168) and MrBayes (v3.2.7a) 
analyses are available as Supplementary Data Files 2 and 3, respectively. Custom R (v3.6.1) code used to plot phylogenetic trees are available 
on github (https://github.com/ThomasAStewart/Qikiqtania). All code is archived at Zenodo (DOI: 10.5281/zenodo.6557684). 


For manuscripts utilizing custom algorithms or software that are central to the research but not yet described in published literature, software must be made available to editors and 
reviewers. We strongly encourage code deposition in a community repository (e.g. GitHub). See the Nature Portfolio guidelines for submitting code & software for further information. 


Data 


Policy information about availability of data 
All manuscripts must include a data availability statement. This statement should provide the following information, where applicable: 


- Accession codes, unique identifiers, or web links for publicly available datasets 
- A description of any restrictions on data availability 


- For clinical datasets or third party data, please ensure that the statement adheres to our policy 


All data are available for download. Computed tomography data sets and STL files of major elements are available for download from MorphoSource (https:// 
www.morphosource.org/projects/000375542). Phylogenetic data are provided in Supplementary Data Files 2 and 3. 
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Ecological, evolutionary & environmental sciences study design 
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Study description The study considers a new fossil specimen, NUFV 137. The specimen was CT scanned to reveal its anatomy. Phylogenetic data, based 
upon published studies and the new CT data, were analyzed to understand the evolutionary relationships of early tetrapods. The 
fossil will be housed to the Canadian Museum of Nature until such time as research and collections facilities are available within 
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Nunavut. 

Research sample This study analyzes specimen NUFV 137. It is the only specimen collected from site NVO401. We describe its collection and its 
morphology. 

Sampling strategy There is only one known specimen (NUFV 137) of Qikiqtania wakei, which we describe in this paper. All fragments have been CT 


scanned, analyzed, and are presented in the study. 


Data collection Computed tomography data were collected at the University of Chicago’s PaleoCT scanning facility. The scanner used was a GE 
Phoenix v|tome|x 240 kv/180 kv scanner. Data were collected by authors TAS and JBL. 


Timing and spatial scale | CT data were collected between Jan. 2020 and Sept. 2022 according to machine availability. 


Data exclusions No data were excluded from analyses. 

Reproducibility There is only one known specimen (NUFV 137) of Qikiqtania wakei; therefore, concerns of replication are not applicable to this 
study. 

Randomization There is only one known specimen (NUFV 137) of Qikiqtania wakei; therefore, concerns of randomization are not applicable to this 
study. 

Blinding There is only one known specimen (NUFV 137) of Qikiqtania wakei; therefore, concerns of blinding are not applicable to this study. 


Did the study involve field work? Xx Yes [| No 


Field work, collection and transport 
Field conditions High Arctic desert with average temps O centigrade with continual wind and occasional snow and sleet. 
Location Bird Fiord, Ellesmere Island Canada. Coordinates of site: N77°10.235’ W86°11.279’. 
Access & import/export — Fieldwork done with permits from Nunavut Research Board; Ministry of Culture, Heritage and Youth, Nunavut; and support from the 


Polar Continental Shelf Program of Canada. Export permit provided by Canadian Museum of Nature, working on behalf of the 
Nunavut Territory. Field permission was granted by the Hamlet of Grise Fiord, Canada; the lviq Hunter and Trappers of Grise Fiord 


Canada. 

Disturbance Disturbances related to campsites include circles of rocks for tents as well as inorganic and organic waste. Surface disturbance at the 
fossil site includes a 8 meter long terrace dug into the rock. The fossil site was infilled, tent rocks moved and wastes removed or 
burned. 


Reporting for specific materials, systems and methods 


We require information from authors about some types of materials, experimental systems and methods used in many studies. Here, indicate whether each material, 
system or method listed is relevant to your study. If you are not sure if a list item applies to your research, read the appropriate section before selecting a response. 


Materials & experimental systems Methods 


n/a | Involved in the study n/a | Involved in the study 
[| Antibodies | ChIP-seq 
[| Eukaryotic cell lines [| Flow cytometry 
[i] Palaeontology and archaeology [| MRI-based neuroimaging 


[| Animals and other organisms 


[| Human research participants 
[| Clinical data 


Dual use research of concern 


Palaeontology and Archaeology 


Specimen provenance The specimen was collected in Canada, Nunavut, southern Ellesmere Island, near the eastern arm of Bird Fiord, Nunavut 


Paleontological Expedition site NVO401, N77°10.235’ W86°11.279’. Export Permit via Canadian Museum of Nature (Ottawa) on 
behalf of Nunavut, August, 2004. 
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Specimen deposition To be returned to the Canadian Museum of Nature upon publication, per agreement of Export and Paleontology permits. 

Dating methods No new dates are provided. 

[| Tick this box to confirm that the raw and calibrated dates are available in the paper or in Supplementary Information. 

Ethics oversight No ethics oversight was required by the Canadian Museum of Nature, Nunavut, or the University of Chicago and Academy of Natural 


Sciences of Drexel University. The field and environmental protocols were approved by the Hamlet of Gris Fiord, Canada; the lviq 
Hunter and Trappers of Grise Fiord Canada, and the Polar Continental Shelf Project of Canada. 


Note that full information on the approval of the study protocol must also be provided in the manuscript. 
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